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ABSTRACT 

In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin 
of the thermal X-ray component may be either release of heat from the core of the 
neutron star, or continuing low-level accretion. In general, heat from the core should 
be stable on timescales < 10^ years, while continuing accretion may produce vari¬ 
ations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, 
Aql X-I) have shown variations in their thermal components on a range of timescales, 
several others, particularly those in globular clusters with no detectable nonthermal 
hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we con¬ 
strain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over 
^ 10 years. We find no evidence of variations in 10 cases, with limits on temperature 
variations below 11% for the 7 qLMXBs without powerlaw components, and limits 
on variations below 20% for 3 other qLMXBs that do show non-thermal emission. 
However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 
CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, 
suggesting the presence of continuing low-level accretion. This work adds to the evi¬ 
dence that the thermal X-ray component in quiescent neutron stars without powerlaw 
components can be explained by heat deposited in the core during outbursts. Finally, 
we also investigate the correlation between hydrogen column density (N//) and optical 
extinction (Ay) using our sample and current models of interstellar X-ray absorption, 
finding A//(cm“^) = (2.81 ± 0.13) x lO^^Ay. 

Key words: accretion, accretion discs, binaries: close, stars: neutron. X-rays: binaries, 
globular clusters: individual: Terzan 5; NGC 6440; NGC 6266 


1 INTRODUCTION 

Transient low-mass X-ray binaries (LMXBs) show occa¬ 
sional outbursts, separated by periods of quiescence (typ¬ 
ically months to years) in which their X-ray emission dra¬ 
matically drops, to Lx ~ 10^^ “10^® erg s“^, and there is lit¬ 
tle or no accretion occurring. The X-ray spectra of quiescent 
LMXBs (qLMXBs) containing a neutron star (NS) consist of 
a soft blackbody-like component (if fit with a blackbody, giv- 
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ing typical temperatures of 0.2-0.3 keV) and sometimes an 
additional harder component, typically fit with a powerlaw, 
with photon-index betwee n 1 and 2 (e.g. ICampana et al.l 
ll998l : lRlitledge et al.lll999h . The origin of the powerlaw com¬ 
ponent in these systems may be continuing low level ac¬ 
cretion onto the NS, an accretion shock between infalling 
matter and the NS’s magnetosphere or a pulsar wind, 
or synchrotron radia t ion from a pulsar wind nebula (e.g. 
Campana et al 1 1 19981: Bogdanov. Grindlav fc van den Berd 
2 OO 5 II . Recent X-rav observations support boundary-layer 
emission from accretion onto the NS as the likely origin of 
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the hard powerlaw component at least in some sources (e.g . 
Cen X-4 and PSR J1023+0038, Chakrabarty et al.| 201-^: 
D’Angelo et al.ir2015l: lArchibald et al.ll2014l:IWiinands et alJ 


20Tj). 


The blackbody-like component can be well-described 
by NS hydrogen atmosphere models, with an implied 
radius consistent with t heore tic ally predicted NS radii (e.g. 
Rajago£al_&_^ommi| Il99fil: IZavlin. Pavlov fc ShibanovI 


19961: iRutledge et alJ Il999l. l2001^ . The soft component 


has been widely interpreted as the slow (10"* years) re¬ 
lease of heat from the core of the NS, deposited through 
deep cru stal heating during previo u s epi sodes of ac¬ 
cretion llSrown, Bildsten &: Rutledgel 1 19981 ). A nearly 
identical soft spectrum can also be created in the up- 
er layers of t he NS at mosphere by low-level ac cretio n 
Zampieri et al.l 1 19951 : IZatie. Turolla fc Treves! I2OO0II , 
though the most detailed model predicts, in addition, an 
optically thin high-temperature bre msstrahlung component 
llPeufel, Dullemond fc Spruid l200ll ~) . whic h may be identi¬ 
fied with the powerlaw component (e.g. ISahramian et aP 
|2014| . However, low-level accretion at sufficiently high 
rates (corresponding to L ~ 10®® erg/s) could maintain 
metals in the atmosphere at abundances sufficient to 
soften the spectrum, leading to an overestimate of R 
ll Rut ledge et al.l l2002bll . The question of the origin of 
the thermal component in qLMXBs is of intense in¬ 
terest, since fitting the thermal spectra of qLMXBs in 
globular clusters is often used to measure the radius of 


tion 

of state (IRutledge et al.l 2002a: Heinke et al. 

l2006a|: 

Webb & Barred 20071: Steiner. Lattimer & Brown 

201ol. 

2013 

: ICuillot. Rutledge & Brown 2011: Servillat et al. 

2012 

: Guillot et al.l l2013l: iGuillot & Rutledgd 20l4 

Lattimer & Steinerl l2014l: Heinke et al.l 20141: Ozel et al. 

2OI5I: see alsolMiller et al.||2013|). Confirmation that thermal 


spectra of qLMXBs (especially those without powerlaw 
components) are not produced by accretion would eliminate 
the possibility that metals remain in the atmosphere, and 
thus eliminate a systematic uncertainty in this method of 
constraining the NS radius. 

At low accretion rates, the NS’s magnetic held is ex¬ 
pected to exercise a propeller effect, w hich retards the accre¬ 
tion of material onto the NS surface (llllarionov fc SunvaevI 
I1975I ). The mechanism is that ionized disc material will be¬ 
come attached to the magnetic field lines when the magne- 
tospheric pressure exceeds the gas pressure; if this bound¬ 
ary (the edge of the magnetosphere) occurs at a radius 
larger than the corotation radius, the material is accel¬ 
erated to velocities higher than the Keplerian orbital ve¬ 
locity, throwing it away, and possibly out of the system. 
On the other hand, simulations of low-level accretion onto 
a magnetic p ropeller predict some material will reach the 
neutron star llRomanova et al. 2()0^ i^lkarni fc Romanoval 


I2OO8I : iD’Angelo fc Spruid I2OIC . l20ld ). Pulsations have 
cently been identified from two (relatively bright) qLMXBs, 
the transitional pulsars PSR J1023-1-0038 at Lx=3 x 10®® 
erg/s, and XSS J12270-4859 at Lx =5 x 10®® erg/s, prov¬ 
ing that matter accretes onto the s urface even though it 
seems to be in the propeller regime llArchibald et al.ll2014l : 
IPapitto et al.l[^14l 'l. Detections of pulsations at low accre¬ 


tion rates from X-ray pulsars (with known magnetic fields) 
in high-mass X-ray binaries also s uggest continued accre¬ 
tion in the propeller regime (e .g. iNegueruela et al.l l2000l : 
iReig. Doroshenko fc Zezaal2014l '). Note that pulsations have 
not yet been dete cted from quiescent L MXBs other than 
transitional MSPs (ID’Angelo et al.l|2015l : Elshamouty et al., 
in prep.). 

Thus, there is evidence that continued accretion 
can produce an X-ray spectrum consistent with the 
soft component in qLMXBs, by heating the N S sur- 


face (Zamgier^jt^ Il995l : IZane. Turolla &: TrevesI 120001 : 
iDeufel. Dullemond fc Spruid |20(5I1 'I and that continued ac¬ 
cretion can continue at very low accretion rates. Testing 
for variability is a promising method to constrain whether 
thermal emission from qLMXBs is driven by accretion. 
X-ray emission driven by heat emerging fro m the core 
shou ld be stable on timescales ^ 10^ y ears (IColpi et al.l 
I2OOII : IWiinands. Degenaar fc Page! l2013l 'l. For NSs that 
have undergone recent episodes of a ccretion, decays of the 


thermal component are seen (e.g. Wiinands^jt^^y_ 200^ 


Cackett et ahll^OSl : iDegenaar et al.ll201ll : Fridriksson et all 


201ll) . attributed to heat leaking from shallower levels i 


the crust, where it was deposited during the ou tburst 
(|Ushomirskv fc Rutledgeir200ll : IRutledge et al.ll20023 L This 
variation should be a monotonic decline on a time scale 
of years/decades. Possibilities for variations of the ther¬ 
mal component without accretion, based on changes in the 
chemical composition of the env e lope, have been advance d 
llBrown. Bildsten fc Chan3 l2002l : IChang fc Bildst^ l2004l l , 
but these do not predict significant changes on timescales 
of years without intervening outbursts. Alternatively, if the 
thermal component is caused by continuous accretion onto 
a NS, the temperature may vary in either direction, on a 
range of timescales. 

Strong evidence for variability of the thermal com¬ 
ponent of some qLMXBs during q uiescence has been re¬ 
ported for the transients Aquila X-]_jRutledgeet_^yj2002^; 


Cackett et al. 


I2OIIIL Cen X-4 (ICampana et al.|[ l997l. |2004 
Cackett et al.l 2010l ). XTE J1701-462 dFridriksson et al l 
2010ll . and MAXI J0556-332 llHoman et al.ll2014l 'l. Each of 


the transient qLMXBs that showed strong variation in the 
thermal component also showed a relatively strong power- 
law component , typically making up ^^50% of the 0.5-10 
keV flux (e.g. Rutledge et al.l l2002bl : ICackett et al.l I 2 OIC 1 I : 
iHoman ef^jal 2014 ). So far, this is consistent with the sug¬ 
gestion ( Heinke et al.l[2003dll that continued accretion is re¬ 
sponsible for b oth the powerlaw component (at least for 
Lx,^10®® erg/s. [jonker et'^l2004| ) and variability. 

Globular clusters (GCs) are highly efficient factories for 
producing X-ray sources (-^lOO times more efficient than 
the r est of the Galaxy, per unit mass; iKatj I1975I : IClarkI 
ll975ll .The densest and most massive globular clusters have 
been shown by Chandra to contain multiple soft X-ray 
sources with Lx in the 10®®-10®® erg/s ran ge, whose spec¬ 
tra strongly indicate that they are qLMXBs (iGrindlav et al 


2001; Poffie^_et_alJl200^ Heinke et al.ll200'^ lM Poolev et al 

2003 ; Maxwell et al.ll2012l) . These can be discriminated from 

other sources (spectrally harder cataclysmic variables, chro- 
mospherically active binaries, and millisecond radio pu l- 
sars) through spectral fitting (e.g. IRutledge et al.l l2002al 'l. 

























































































































































































Thermal Lx variations in qLMXBs 3 


or through X-ray colours and luminosities l|Heinke et alJ 
l2003dl : fPoolev &: Hutll2do^ '). The majority of the candidate 
qLMXB s show no evidence fo r powerlaw components in their 
spectra (|Heinke et alJlioO.ldh . This may be a selection effect, 
since it is easier to discriminate qLMXBs from other sources 
if they show purely thermal spectra. Deep Chandra obser¬ 
vations of globular clusters identify a significant population 
of candidate qLMXBs with strong powerlaw components, 
which require high-quality X- ray spectra to confidently iden¬ 
tify t he thermal components (iHeinke. Grindlav &: Edmonds! 
I 2 OO 5 I I. and half of the transient cluster LMXBs followed up 
with Chandra are too faint, or too dominated by a pow- 
erla w component, for detection of a thermal component 
(e.g.lWiinainds et al.ll2005l : iHeinke et al.ll201ol : iLinares et al.l 

I2OI4IL 


Several qLMXBs in globular clusters with thermal 
components have been examined to search for varia¬ 
tions between observations. Several show flux variations, 
but spectral analyses have permitted these variations to 
be driven only by ch anges in the normalization of the 
powerlaw componen t (|Heinke ^_GrjrKllja,v fcEdmondsI l2005l : 
ICackett et al.l I 2 OO 5 I : iBahramian et ahT 2014 1. or of ob- 
scuri ng material-the latter especially in edge-on systems 
fe.g. iNo wak. Heinz fc BeeelmarJ 1200IHeinke et al.lf^OSd : 
IWiinands et al. 2003ll . A few qLMXBs without detectable 
powerlaw components have repeated deep observations with 
Chandra, permitting sensitive measurements for variabil¬ 
ity; these include (with 90% confidence va riability lim¬ 
its) 47 Tuc X7, AT/T <1.0% over 2 years llHeinke et al.l 
l2006al): NGC 6397 U18, Ar/r <1.4% over 10 years 
i Guillot. Rutledge fc Brownll201ll: Heinke et al. 201^; M28 


source 26 (no variation found. IServillat_et_aJjj201^ and uj 
Cen, ATjT <2.1% over 12 years ( Heinke et al. 20141) . 


We are interested in whether this lack of variabil¬ 
ity is the norm in thermal components of globular cluster 
qLMXBs, and whether (if it is observed anywhere) there 
is a correlation with the presence of powerlaw components. 
Strong evidence that globular cluster qLMXBs without pow¬ 
erlaw components are not variable would increase conhdence 
in the assumptions used to derive radius measurements for 
their NSs and obtain constraints on the dense matter equa¬ 
tion of state. We choose to focus on qLMXBs in globular 
clusters, since globular clusters are the only places where we 
can find and study large populations of qLMXBs in single 
Chandra pointings (and with relatively low optical extinc¬ 
tion, compared to many qLMXBs in the Galactic Plane) and 
distances to qLMXBs in globular clusters are known better 
than to the ones in the field. 


In this paper, we analyze a sample of 12 qLMXBs in 3 
globular clusters, NGC 6266 (M 62), NGC 6440, and Terzan 
5, searching for variations over multiple observations span¬ 
ning roughly a decade. Each of these dense globular clus¬ 
ters have at least 4 candidate qLMXBs, and multiple C/ian- 
dra/A ClS observations. When this work w as in an advanced 
state, IWalsh. Cackett fc Bernardinil (l2015l ) reported on the 
analysis of 9 similar sources in Terzan 5 and NGC 6440 as 
considered here, though using fewer Terzan 5 observations, 
and not including NGC 6266. We directly compare our re¬ 
sults with those of Walsh et al., finding general agreement. 
In U21 we describe our data reduction and analysis methods. 


In 331 we show the results of our analysis, and in 311 we 
discuss the implications. 


2 DATA REDUCTION AND ANALYSIS 


All datasets were obtained using ACIS-S in Faint, Timed 
Exposure mode. The data covers a time span of 9 to 12 
year s for the GCs in our s ample (Table [J). We used CIAO 
4.6 (iFruscione et al.l l2006l) and CALDB 4.6.2 for data re¬ 
processing following standard CIAO science thread^ We 
chose known (and candidate) LMXBs with more than ~ 60 
photons (per epoch) for spectral analysis, resulting in a total 
sample of 12 sources. These targets are tabulated in Tabled 
We performed spectral extraction using the task specextract 
and performed spectral analysis in the 0.3-10 ke V energy 
range using HEASOFT 6.16 and XSPEC 12.8.2 jArnaud 
Il996l) . We combined (using the CIAO dmmerge tool) spec¬ 
tra from observations which occurred within a month for 
our targets in Terzan 5, to improve the spectral quality. We 
grouped each spectrum by 15 counts per bin, and used chi- 
squared statistics for analysis. Details of the spectral extrac¬ 
tions and analyses for each cluster are discussed in 33 All 
uncertainties reported in this paper are 90% confidence. 

We fit the spectra with an absorbed neutron star at- 
mosp here, plus a powerlaw (if needed). We used the NSAT- 
MOS jHeiiike et al.l[^06a^ NS atmosphere model, with the 
NS mass and radius frozen to the canonical values of 1.4 Mq 
and 10 km, the normalization fixed to 1 (implying radiation 
from the entire surface), and the distance fixed to the esti¬ 
mated distance of the cluster (see 33- In cases where this 
single-component model does not produce a good fit, leav¬ 
ing significant residuals at high energies(^3 keV), a pow¬ 
erlaw component (PEGPWRLW, with normalization set to 
represent flux in the 0.5-10 keV range) was added to the 
model. In cases where it was unclear whether the fit was 
improved by adding a second component, we checked this 
by an F-test. We set the powerlaw photon index (P) as a 
free parameter when fitting datasets where the data qual¬ 
ity above 3 keV is high enough to constrain this quantity. 
However, in datasets where the data quality is not sufficient 
to do this, we fixed the photon index to 1.5, which is a 
typic a l value from LMXBs in quie s cence dCampana et al.l 


1T9^ 2004 Fridrikssori et all l201ll : IPeaenaar et al. 2011 
Chakrabartv et al. \\2Qli)- 


We fit all spectra available for each source simultane¬ 
ously. We tied the N_tf and powerlaw photon index to a 
single value for each source, and let the NS temperature 
and powerlaw flux vary between observations. There is good 
evidence that Ah does not vary during the spectral evo¬ 
lution of most LMXBs llMiller. Cackett fc Reidl2009l ). Al¬ 
though correlated variations in N/r and powerlaw photon 
index have been proposed to explain spectral var iability 
in Aql X-1 in quiescence (ICamoana fc Stellal 12003 ). such 
a model requires the powerlaw to dominate over the ther¬ 
mal component of the spectrum, and forces the photon in¬ 
dex to (unreasonably) high values (^ 2.5). Fitting individ - 
ual high-quality spectra of Aql X-1 (ICackett et al.ll201ll ). 


^ http://cxc.harvard.edu/ciao/threads/index.html 
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SAX ,71808.4-3658 dHeinke et alJl2009l l. and Terzan 5 X-3 
dBahramian et ani2014l 7 in quiescence gives consistent val¬ 
ues (within the errors) between epochs for both N^f and 
the powerlaw photon index. Deep observations of Cen X- 
4 do show ev i dence for powerlaw photon index changes 
dCackett et ahlBoid l. changing from ~1.7 in all observa- 
tions where Lx < 2 x 10^^ erg/s, to -^1.4 at h i gher L x values 
dChakrabartv et al.ll2014l 7 . iFridriksson et akl (120101 7 also find 
evidence for variation in the powerlaw photon index in XTE 
J1701-462, but only when Lx reaches > 2 x 10^^ erg/s, well 
above the range discussed here. Summarizing these observa¬ 
tions, when the total Lx does not vary by a large factor, the 
powerlaw photon index appears to remain roughly constant. 
As we will see below, our observations do not reveal large 
swings in Lx, so we feel that keeping a fixed powerlaw index 
is justihed. 

As we are looking for relative variations in kT, to es¬ 
timate the uncertainties in the NS temperature and pow¬ 
erlaw flux, after initially fitting the Nh we then froze N_tf 
to its best-fit value. In cases where the source shows signif¬ 
icantly higher N_tf value than the average of the cluster, we 
investigated the possibility of intrinsic absorption leading to 
variations in Njf (!l3]). 

We also search for signs of temporal variations during 
each observation for all sources. (Due to the limited sig¬ 
nal, we explore only total flux variations.) We extracted 
background-subtracted lightcurves in 0.3-7 keV band for 
each source using CIAO task dmextract. Due to the faint 
nature of the sources in our study, we binned each lightcurve 
by 1000 s. We used the FTools task Icstats to perform chi- 
squared and Kolmogorov-Smirnov (KS) tests for variations. 
Since we ran these tests on 80 lightcurves (all sources in our 
sample in all observations useqj), we require a false-alarm 
probability below 5%/80=0.06% to identify variability at 
95% confidence (Note that null-results from a test do not 
prove the source did not vary). 


3 RESULTS 

3.1 Hydrogen column density 


We constrained the hydrogen column density {Nh) for each 
GC using our sample of LMXBs. In our primary analysis in 
this paper, we forced the N_tf to a single value for all obser¬ 
vations of each source. We used th e Tuebingen-Boulder ISM 


absorp tion model (TBA BS) with Wihng^^Allen__fc_JvlcCray 


(200C) abundances and iBalucinska-Church fc McCammon 
(1992) cross sections. fWe tested the IVerner et al.l ( 199^ 
cross sections, and found no diff erence in our results). We 
also tried using the abundances of lAnders fc Grevessd (Il989l 7 
for comparison (see below). After constraining 'Nh values for 
each target, we calculated the mean value for sources in each 
cluster. (We excluded Terzan 5 GX 9 from this analysis due 
to possible intrinsic absorption-see il3.4l 7 We compare these 
values with those produced by adopting reddening values, 


^ We excluded one of Terzan 5 observations (Obs.ID 13705) from 
this investigation due to its short exposure (14 ks). 


E(B-V), from the Harris catalog (lHarrislll99d 2010 revi¬ 
sion) and the relation between Nh and optical extinction 
(Av) from lGiiver fc 6z'dll2009l 7: 

Nnicm-^) = (2.21 ± 0.09) x lO'^^Av (1) 


(Table [3]). The Nh value s calcu lated from Ay using the 
relation of iGuver &: Ozell (l2009l7 agree nicely with fits to 
the X-ray spectra using the lAnders fc Grevessel (19^) 
abundances-which makes sense, as iGtiver fc Ozell (2009) 
used Nh estimates from the literature on X-ray studies of 
supern ova remnants, most of which used lAnders fc Grevessel 
(Il989l7 abundances. How ever, Ng^ values measur ed from 
spectral htting using IWilms. Allen fc McGravl (l2000l) 
abundances a r e ^2 7 % higher than predicted using 
Gfiver fc Ozel (l2009l). This makes sense, because the 
Wilms. Allen fc McGravl (l2000l ) abundances are typically 
~30% lower for X-ray relevant elements. We therefore intro¬ 
duce a new relation between Nh and Ay, des i gned specifi- 
cally for use with the I Wilms. Allen fc McGravl I 2 OOOI ) abun¬ 
dances. 


Nnicm-'^) = (2.81 ± 0.13) x lO^Mv 


( 2 ) 


where the error is the 1-a scatter we measure in the rela¬ 
tion from the fit. We derived this equation from fitting a 
linear model to Nh values measured for 11 out of 12 sources 
in our sample (we excluded Terzan 5 GX 9 due to possi¬ 
ble presence of enhanced absorption). We also investigated 
the effects of including du st scattering for poin t sources by 
adding a scattering model (IPredehl et al.|[2003l) to our spec¬ 
tral model, with Ay fixed to 0.117xAfir/(10^^ cm“^). We 
found that dust scattering causes a slight decrease of «!% 
in the Nh values, and thus increases the Ay-to-Nn con¬ 
version f actor to 2.84x10^^ cm~^. As this work was being 
Hnalized, iFoight et al.l ([201^ published their investigation 
of the Ay-to-Nn correlation using optical vs. X-ray studies 
of sup ernova remnants, also using IWilms. Allen fc McGravl 
(l2000l7 abundances, which found the conversion factor to be 
2.87T0.12 X 10^^ cm“^. This is in complete agreement with 
our results. 


3.2 NGC 6266 

NGC 6266 (M 62) is a massive dense globular cluster. At 
least 26 X-ray sources can be identihed in Chandra images 
of the cluster, of which 5 have been suggested to be qLMXBs 
on the basis of their soft spectra and X-ray luminosities 
dPoolev et al.ll200^ : iHeinke et al.ll2003(j7 . One of these (the 
brightest, and with the hardest X-ray spectrum; reasonably 
ht by an absorbed powerlaw of photon index 2.5T0.1) shows 
strong evidence (f rom its radio/X-ray flu x ratio) for contain¬ 
ing a black hole (iGhomiuk et al.ll2013l 7. We use two deep 
Chandra/ACIS observations (Table [TJ to study the other 4 
soft candidate qLMXBs, which we name GX 4, GX 5, GX 
6 , and GX 16 (ordering the detected sources by Chandra 
countrates). In the 2002 observation, the cluster core is 2.5 
arcminutes off-axis, and the sources have visibly elongated 
point-spread functions. To address this issue, we used ellip¬ 
tical extraction regions of the same size for all sources. The 
2014 observation was performed on-axis, so we used circular 
extraction regions. We show images of these observations 
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Target 

Obs. ID 

Date 

Exposure (ks) 

NGC 6266 

02677 

2002-05-12 

62 


15761 

2014-05-05 

82 

NGC 6440 

00947 

2000-07-04 

23 


03799 

2003-06-27 

24 


10060 

2009-07-28 

49 

Terzan 5 

03798 

2003-07-13 

40 


10059 

2009-07-15 

36 


13225 

2011-02-17 

30 


13252 

2011-04-29 

40 


13705 

2011-09-05 

I4A 


14339 

2011-09-08 

CO 


13706 

2012-05-13 

46 


14475 

2012-09-17 

30 


14476 

2012-10-28 

29 


14477 

2013-02-05 

29® 


14625 

2013-02-22 

49® 


15615 

2013-02-23 

84® 


14478 

2013-07-16 

29 


14479 

2014-07-15 

29*^ 


16638 

2014-07-17 

72*^ 


15750 

2014-07-20 

23*^ 


Table 1. Chandra ACIS observations used in this study. A, B, C: We merged spectra extracted from all observations marked A (and 
merged those marked B, etc.), as they occurred close in time. 


Source 

CXOGlb J 

Lx 



(lO^^erg s“^, 0.5-10 keV) 

NGC 6266 CX4 

170113.09-300655.43 

4.5±0.3 

NGC 6266 CX5 

170113.10-300642.33 

4.2±0.3 

NGC 6266 CX6 

170113.76-300632.48 

3.9±0.3 

NGC 6266 CX16 

170112.58-300622.38 

1.0±0.2 

NGC 6440 CXI 

- 

10±1 

NGC 6440 CX2 

- 

14±1 

NGC 6440 CX3 

- 

8±1 

NGC 6440 CX5 

- 

7.4±0.8 

Terzan 5 CX9 

174804.8-244644 

9.0±0.9 

Terzan 5 CX12 

174806.2-244642 

7.0±0.8 

Terzan 5 CX18 

174805.2-244651 

6.3±0.6 

Terzan 5 CX21 

174804.2-244625 

4.4±0.5 


Table 2. Candidate LMXBs studied in this paper. Luminosities are based on latest epoch used for each GC and are measured in this 
study. Reported uncertainties for lum inosities are 90% c onfidence and d o not include uncert ainties in distance. References for source 
identification: This work (NGC 62661. [Poolev et al.ll200^ (NGC 6440), fOl lHeinke et al]|2006bl (Terzan 5). 


and extraction regions in Fig. [T] We ran wavdetect to de¬ 
termine accurate positions and brightnesses of these X-ray 
sources, and corrected the astrometry by applying coordi¬ 
nate shifts to match the radio and X-ray positions of M62 
CXI. 

We assumed a distance of 6.8 kpc, and a redden¬ 
ing of E(B-V)= 0.47 mag, for NGC 6266 in our analy¬ 
sis dHarri 3 Il996l . 2010 revision). Assuming Ry=3.1, this 
reddening giv e s Av =1.46 mag. The Ay / Nh relation of 
ICuyer fc Ozell ll2009h which is based on I Anders fc Grevessd 
( ^98^ abundan ces, predicts Njf= 3.2x10^^ cm ^ for this 
cluster based on lAnders fc Grevessel (Il989l ) abundances; in 


i|3.1l we use N_fr values measured from our spectral analysis 
of X-ray observations to derive a new relationship between 
Nh and Ay, Nh= 2.8x10^^ Ay , specifically for use with the 
IWilms. Allen fc McGravl (l2000l) abundances. 

Due to the relatively low interstellar absorption and rea¬ 
sonable fluxes of our targets, the spectra from this cluster 
are of high quality (Fig. [2] top left). We found no evidence of 
temperature variations in these objects (Tabled. For three 
of our targets (CX4, GX5 & CX6) we constrain temperature 
variations to less than 5%. For CX16 this limit is < 10% due 
to its faintness. 

NGG 6266 GX4 is the only candidate qLMXB in NGG 
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GC 

E(B-V) 

(mag) 

Av 

(mag) 

Calculated Nh (Guver09) 
(xlO^^ cm“^) 

Measured N^f (Anders89) 
(xlO^^ cm~^) 

Measured Nh (WilmsOO) 
(xlO^^ cm“^) 

Discrepancy 

(%) 

NGC 6266 

0.47 

1.46 

3.2±0.2 

3.2±0.4 

4.1±0.4 

22±9 

NGC 6440 

1.07 

3.32 

7.3±0.5 

7.2±0.6 

10±1 

27±9 

Terzan 5 

2.61 

8.09 

18±1 

17±2 

26±2 

31±9 


Table 3. Comp aring; Nh for GCs ba sed on different assumptions. Calculated Nh is calculated base d on relation between and 
provided by iGiiver &; Ozdl ||2009|) . The uncertainties reported for calculated values based on iGiiver &: Qzdl ll2009h are only 
statistical uncertainties from the correlation. Measured Njf values are based on X-ray spectroscopy of the sources in our sample (averap^e 
and u nce rtainty calculated for each GC ; see text) with two different assumptions for the abu ndances of e lements jAnders Sz Grevessd 
Il989l and lWilrns. Allen &: McCravll2000h . E(B—V) values are obtained f rom the Harris Catalog llHarris|[l996t) : to obtain Ay, we assumed 
Rl/=3.1. Discrepancy is calcu lated between the m easured Njy using IWilms, Allen &: McCravl ll200Cli) abundances and the calculated 
values based on the relation of lGitver &: OzeH All uncertainties are 90% confidence. 


6266 that shows evidence for a powerlaw component (Ta¬ 
ble 14)1. To constrain the contribution from the powerlaw, we 
fixed the photon index to 1.5 due to the limited spectral 
quality at higher energies. We notice that there is a hint 
of a faint source 1" from CX4 (southwest of CX4, vis¬ 
ible when the image is over-binned). To address this, we 
performed a second set of spectral analyses with a smaller 
source extraction region excluding the region of the possible 
nearby source, but we found no significant spectral changes 
suggesting the non-thermal emission is therefore intrinsic to 
CX4 (Table |4] contains fits performed with the standard ex¬ 
traction regions). 

We searched for short (~1000 s) timescale temporal 
variations in each observation by performing chi-squared and 
KS tests, and found no evidence of variation for any of our 
targets in NGC 6266. No source shows probability of con¬ 
stancy less than 0.06% which is the 95% confidence range 
when considering all trials (80, for all sources in this study). 


3.3 NGC 6440 


NGC 6440 is a moderately extinguished GC near the 
Galactic center. There are at least 25 X-ray sources 
identified in this cluster, of which two are known tran- 
sient LMXBs (NGC 6440 CX 1 = SAX ^748.9-2021 
indZandet^ 199^ indZandetay 200 ll : iPoolev et al'l 


2 OO 2 I: IBozzo, Kuulkers fc Ferrign3 l2015l: NGC 6440 X-2, 
Altamirano et al. I I 2 OIC 1 I : iHeinke et al. [ioi^. This cluster is 


8.5 k pc away and has a reddening of E(B-V) = 1.07 (iHarrisI 
Il996ll . We used three long Chandra/ACIS observations of 
this cluster, obtained in 2000, 2003, and 2009. We study the 
four brightest identified qLMXB candidates in or near the 
core of this cluster, CX 1, CX 2 (not to be confused with 
the AMXP denoted as NGC 6440 X-2 , whic h is not detected 
in deep quiescence, e.g. iHeinke et al1l2010l l. C X 3, and CX 
5, wh ich were identified as likely qLMXBs bv IPoolev et al.l 

ll2002ll . 


The Chandra/ACIS observation of NGC 6440 in 2009 
occur red while NGC 6440 X-2 was in outburst llHeinke et al.l 
l2010lf and the image is severely piled-up in the vicinity of 
this source. However, we were nevertheless able to use this 
observation for spectral analysis of qLMXBs, as NGC 6440 
X-2 is located further away from the core of the cluster (~ 
12 ”) and the outburst was relatively faint (peaking at ~ 
1 X 10^® erg s“^). We addressed the contamination from the 


X-2 outburst on affected sources (CX 1, CX 2 & CX 3) by 
extracting background spectra from multiple regions at the 
same angular distance as the target from X-2. Our targets 
and chosen extraction regions are shown in Fig. [3] 

Two of the sources in this cluster (CX 1 and CX 3) 
show excesses at high energies to the NSATMOS fits, so we 
added a powerlaw component to their spectral models. In 
both cases, we froze the powerlaw photon index to the value 
found in the best fit to all three observations. We also note 
that there is a hint for the presence of a very faint source 
near CX 3 (~10 times fainter than CX 3, south of it, visible 
when the image is over-binned), however it does not affect 
the spectrum. The results of our spectral analyses for NGC 
6440 a re tabulated in Table 1^ 1 Walsh. Cackett fc Bernardinil 
ll2015ll measure powerlaw contributions consistent with our 
results for CX 2 and CX 3 in both 2003 and 2009 observa¬ 
tions; We note that our extraction regions might differ fr om 
the ones chosen bv lWalsh. Cackett fc Bernardinil (l2015l l as 
they do not discuss their choice of extraction regions. This 
is specially important due to the large and spatially varying 
background from the outburst of NGC 6440 X-2. 

From our sample, CX 2 and CX 3 show no evidence of 
variation in the NS t emperature. This is in genera l accor d 
with the results of IWalsh. Cackett fc Bernardinil ll2015h . 
These authors include a short (4 ks) Chandra/ACIS obser¬ 
vation of NGC 6440, 12 days after our 3rd observation, and 
find that the flux increases in this observation for CX 5. If 
the spectrum is assumed to consist of only a NS atmosphere 
model, then the NS temperature must have incre ased in CX 
5. However. IWalsh, Cackett fc Bernardinil (l2015h note that, 
since CX 5’s spectrum may contain a powerlaw in this ob¬ 
servation, the variability cannot be conclusively attributed 
to changes in the NS temperature. We find no evidence for 
variation in the spectrum of CX 5 between 2000 and 2009. 
We also detect the presence of a weak powerlaw component 
in the 2000 epoch at 1.7 a. 

We found possible evidence for thermal variations in 
CX 1 in quiescence (Fig. (2] top right). We significantly de¬ 
tect a powerlaw with photon index 1.7±2 in the 2000 ob¬ 
servation, but not in the 2003 or 2009 observations. Fixing 
this powerlaw index to the best-fit index, our best fit finds 
a — 10lg% change in temperature from 2000 to 2003, with 
the 2009 observation consistent with the 2003 observation 
(Table [ 5 ]). This marginal evidence for variation suggests the 
presence of continuing accretion during the 2000 observa- 
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Figure 1. X-ray (0.3-7 keV) images of NGC 6266 as seen by Chandra-ACIS in 2002 (left) and 2014 (right). Source extraction regions are 
represented by solid circles/ellipses, and background extraction regions are represented by dashed circles. In the 2002 observation, the 
centre of NGC 6266 was located 2.5 arcminutes off-axis. To account for the distortion of the point-spread function this offset induced, 
we used elliptical extraction regions for that observation. 


Source 

Nh (lO^l cm-2) 

Year 

logT (K) 

% Variation 

% Powerlaw flux fraction (0.5-10 keV) 

NGC 6266 CX 4 

3.9±0.4 

2002 

6.057±0.007 

- 

8±7 

Powerlaw r=(1.5) 
xl/d.o.i= 0.87/56 


2014 

6.056±0.006 

< 3.2 

12±5 

NGC 6266 CX 5 

3.9±0.4 

2002 

6.055±0.007 

— 

< 9 

Powerlaw r=(1.5) 
xl/d.od= 0.87/52 


2014 

6.062±0.006 

< 4.7 

< 7 

NGC 6266 CX 6 

4.2±0.4 

2002 

6.056±0.006 

_ 

< 9 

Powerlaw r=(1.5) 
xl/d.od= 0.64/49 


2014 

O.UOO_Q QQ.J, 

< 3.6 

< 7 

NGC 6266 CX 16 

4.3±0.6 

2002 

5.96±0.01 

— 

< 23 

Powerlaw r=(1.5) 
xS/d.o.f = 1.21/14 


2014 

5.94±0.01 

< 9 

< 15 


Table 4. NGC 6266: Results of spectral analysis using an absorbed NS atmosphere (TBABS*NSATMOS), plus a powerlaw added when 
needed. All uncertainties are 90% confidence. Temperature variations are calculated from the first observation in 2002. Powerlaw fractions 
are fractions of the total 0.5-10 keV flux. 


tion (plausible considering that CX 1 erupted into outburst 
two years before, and one year after, the 2000 observation, 
lin’t Zand et al.ll200ll ). 

ICackett et al.l (l2005l ) investigated the spectral varia¬ 
tions of CX 1 using the first two Chandra/ACIS obser¬ 
vations of NGC 6440, but they argued that the observed 
spectral variations were likely due to changes in the pow¬ 
erlaw component exclusively. Cackett et al’s best fit uses a 
powerlaw component with photon index 2.5(±1.0), which 
explains the difference with our best fit (photon index 
1.7±2); a softer powerlaw can explain flux differences in 
the soft X-rays, whereas our harder value of photon in- 
dex require changes in the thermal component to explain. 
IWalsh, Cackett fc Bernardinil ll2015h also fit the spectra of 
CXI in these 3 observations (adding the short 2009 obser¬ 


vation mentioned above) using C-stats in Xspec, find a dif¬ 
ferent value from I Cackett et al.l (l2005l) and basically agree 
with our results-the powerlaw photon index is measured at 
1.6±1.0, and the 2000 data indicate a slightly higher NS 
temperature compared to the 2003 data, at la confidence 
(though not at 99% confidence). Walsh et al. also conclude 
(from fitting) that it is not clear whether the variation in 
CXI was due exclusively to changes in the powerlaw com¬ 
ponent, exclusively to changes in the thermal component, 
or to changes in both. We agree with Walsh et al. that it 
is not clear whether the thermal component varied between 
the 2000 observation and later observations. 

Searching for short (~ 1000s) timescale variations in 
each observation for all sources, we found no evidence of 
variation. No source shows probability of constancy less than 
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0.5 1 2 5 10 

Energy (keV) 



0.5 1 2 5 10 



0.5 1 2 5 10 

Energy (keV) 


Figure 2. Examples of extracted spectra; top panels show data (crosses) and model (line). Different colours represent different epochs. 
Note that in addition to changes in the source, changes in the detector may alter the observed spectrum (these changes are included in 
the model). Top left: NGC 6266 CX4 (see ^3.211 . Only two Chandra/ACIS epochs are available, but the spectral quality is high. Top 
right: NGC 6440 CXI (see H3.3D . The black (2000) spectrum is brighter at all energies, particularly at higher energies; we attribute this 
to a difference in the powerlaw component, and possibly the thermal component. Bottom left Sz right: Terzan 5 CX9 Sz Terzan 5 CX12 
respectively, including spectra from 11 epochs (see The spectral fit for CX9 shown here uses a hydrogen atmosphere model. Note 

that there are still clear waves in the residuals (down at 2 keV, up at 3, down again at 5 keV), indicating problems with the fit (though 
the is formally acceptable, 87.42 for 87 dof). These residual patterns are not completely eliminated in a helium atmosphere model, 
though they are reduced. CX12 shows evidence for a strong powerlaw component in this fit, as well as evidence for variation in soft 
X-rays (particularly in the 4th observation, taken in April 2011, colored dark blue here). 


0.06% which is the 95% confidence range when considering 
all trials (80, for all sources in this stndy). 


3.4 Terzan 5 


Terzan 5 is a massive, highly extincted GC near the Galac- 
tic center. It is located at a distance of 5.9±0.5 kpc 
dValenti, Ferraro fc Ong Un|[2007 ) and it h arbours more than 
40 X-ray sources ( Heinke et al.l l2006bll . There are three 
known transi ent LMXBs in this clu ster, the most in any 
Galactic GG dBahramian et aDl2014h . Terzan 5 has shown 
numerous outbursts (an ov erview up to 2012 is given in 
IPegenaar fc Wiinandsl2012l ) , several of which were not accu¬ 
rately localized (the most recent detected outburst without 
accurate localization being in 2002). It is therefore quite pos¬ 
sible that other sources have been in outburst in the recent 
past. 

We use all available Chandra/ACIS observations of this 
GG in which all sources are quiescent to study the brighter 


thermally-dominated qLMXBs without observed outbursts. 
We start wi t h the sample of candidate qLMXBs from 
I Heinke et al.l (l2006al ). which identified thermally-dominated 
qLMXBs by spectral fitting, or by hardness ratio plus 
inferred X-ray luminosity. Analyses of the quiescent be¬ 
haviour of the outbursting sources Ter 5 X-1 (EXO 1745- 
248), X-2 (IGR J17480-2446) and X-3 (Swift J174805.3- 


244637) have been reported in 

Deeenaar Sz WiinandsI 

2011 ; 

Degenaar. Brown & WiinandsI 

20111: Degenaar et al. 

2013; 

Bahramian et al.ll2014l: iDegenaar et al.ll2015l. so we exclude 


those sources here. 


For the purpose of imaging and source detection, we 
merge all available observations following CIAO thread^ 
We corrected relative astrometry in all frames by matching 
the coordinates of 6 bright sources. We ran reprojecCaspect 
to create a new reprojected aspect solution for each obser¬ 
vation. Then we applied these aspect solutions to the event 

® http://cxc.harvard.edu/ciao/threads/combine/ 
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Figure 3. NGC 6440 as seen by Chandra-ACIS in the 0.3-7 keV band in 2003 (left) and 2009 (right). Solid circles represent the source 
extraction regions and dashed circles show the se lected background regions. In the 2009 observation, the image is contaminated by an 
outburst from NGC 6440 X-2 llHeinke et ahll^Oinh . We address this contamination for affected sources by choosing multiple background 
extraction regions at the same distance from X-2 as the target. 


Source 

Nif (10^2 cm-2) 

Year 

logT (K) 

% Variation 

% Powerlaw flux fraction (0.5-10 keV) 

NGC 6440 CXI 

l.liO.l 

2000 

6.16±0.01 


26±10 

Powerlaw r=(1.7) 


2003 

O--l-^_0.02 


< 10 

x5/d.o.f = 0.84/37 


2009 

6.14±0.01 

< 11 

< 8 

NGC 6440 CX2 

l.liO.l 

2000 

6.16±0.01 

_ 

< 6 

Powerlaw r=(1.5) 


2003 

6.179±0.009 

< 9 

< 9 

X?/d.o.f = 0.87/47 


2009 

6.172±0.007 

< 7 

< 7 

NGC 6440 CX3 

0.9±0.1 

2000 


- 

35±17 

Powerlaw r=(2) 


2003 

6.10±0.02 

< 17 

33±14 

xS/d.o.f = 0.57/35 


2009 

6.10±0.02 

< 17 

19±15 

NGC 6440 CX5 

l.OiO.l 

2000 

6.09±0.01 

— 

18±15 

Powerlaw r=(1.5) 


2003 


< 10 

< 10 

xS/d.o.f = 1.00/25 


2009 

6.11±0.01 

< 10 

< 7 


Table 5. Spectral analysis of qLMXBs in NGC 6440. We fit the spectra with an absorbed NS atmosphere (TBABS*NSATMOS), and 
included a powerlaw component when necessary. All uncertainties are 90% confidence. Temperature variations are calculated based on 
the coldest measured temperature for each source. Powerlaw fractions are fractions of the total 0.5-10 keV flux. In cases where the 
spectral quality allows us to constrain the powerlaw photon index (CXI and CX3), we tied the index between observations and searched 
for the best-fit value, then froze the index to that best-fit value; otherwise we fixed the index to 1.5. The photon index used is listed in 
column 1. 


files by running reproject-events, and finally we combined 
the reprojected event files by running reproject-obs. 

We selected the candidate qLMXBs which had sufficient 
counts for spectral analysis in each epoch 60; i.e. more 
than 4 spectral bins of 15 counts), leaving us with a sample 
of four thermally-dominate d qLMXBs, identifie d as CX9, 
CX12, CX18, and CX21 in iHeinke et al.l (l2006al ). We used 
the merged image to choose source and background spectral 
extraction regions. This is important since in cases like CX9 
and CX18, the vicinity of the source is complex and there 
can be multiple sources of contamination. In cases where 
the vicinity of the source is complex/crowded, we extracted 
background spectra from annuli around the targets, exclud¬ 


ing detected sources. The combined image of Terzan 5, with 
our extraction regions, is presented in Fig. |T] 

We combined spectra from observations which occurred 
less than a month apart (using comhinespectra), to obtain 
better statistics. Investigating these observations individn- 
ally, no spectral variation was detected. These observations 
are marked in Table [T] All our targets in Terzan 5 show a 
high energy excess suggesting the presence of a powerlaw 
component in their spectra with average fraction of ~10% 
of the total flux for CX9, CX18, and CX21 and an aver¬ 
age fraction of ~33% for CX 12. Comparing these with our 
targets in NGC 6266 and NGC 6440, where the powerlaw 
component is mostly detected when flux fraction is ^10%, 
suggests that our detection of a powerlaw component in all 
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our targets in Terzan 5 might be due to the greater depth 
of these observations, combined with our relatively higher 
sensitivity to hard vs. soft photons in this heavily absorbed 
cluster. We froze the powerlaw photon index to 1.5 for CX9, 
CX18 and CX21, due to a lack of sufficient signal to con¬ 
strain it. For CX12 we freeze it to the value of 1.8, obtained 
from the best fit to all observations. 

The values of N_tf for CX12, CX18 and CX21 are all 
in agreement with the value of 2.6±0.1 x 10^^ cm“^ from 
previous studies (see 3330. However, CX9 has a hydro¬ 
gen column density of 3.1±0.2 x 10^^ cm“^, higher than 
the measured value for the cluster. We investigate the high- 
resolution reddening map of Terzan 5 (iMassari et al.ll2012l . 
available online from Cosmic-Lab 0 ), and find no significant 
difference (<1.5%) in reddening between the direction of 
CX9 and the rest of the sources in our sample. This suggests 
that the inferred high Njr value may be caused by intrinsic 
absorption in the system. Spectra of CX9 are shown in Fig. 
[3] (bottom left). 

We also performed spectral fits on CX9 by replacing 
the non-magnetic hydrogen atmosphere model (NSATMOS) 
with a non-magnetic heli um atmosphere model (NSX in 
Xspec, fc Heinl^ I 2 OO 9 I L This fit gives a slightly better 
compared to the hydrogen atmosphere model (x^ de¬ 
creased from 114 to 109, where both cases have 102 de¬ 
grees of freedom), and also a lower absorption value of 
Nh = 2.9 ± 0.2 X 10^^ cm“^, consistent with the rest of 
the cluster (Table [S|. These results indicate that a helium 
atmosphere is a possible explanation for the unusual aspects 
of CX9’s spectrum (see Note that there are still clear 
waves in the residuals (down at 2 keV, up at 3, down again 
at 5 keV), indicating problems with the fit (though the x^ 
is formally acceptable, 87.42 for 87 dof). These residual pat¬ 
terns are not completely eliminated in a helium atmosphere 
model, though they are reduced. 

We found no evidence regarding spectral variations be¬ 
tween observations in CX 9, CX 18 and CX 21. However CX 
12 shows marginal evidence of variation over time (Table [T]). 
Our main fit, in which only the NS temperature and pow¬ 
erlaw flux are free among observations, suggests that these 
variations are caused by changes in the temperature of the 
NS, particularly a decrease in the 2011-4 epoch, and to a 
lesser degree in the two 2013 observations. However, the sub¬ 
stantial powerlaw component, which is particularly strong in 
the epochs with the coldest NS temperature measurements, 
suggests that subtle changes in the powerlaw (both flux and 
index) might cause the variations (Fig. [21 bottom right). We 
note that thermal variations are only suggested in 2 epochs 
(2011-4 and 2013-2) out of 11. However we get an acceptable 
chi-squared when fitti ng the temperatures to a straight line . 
Thus we agree with IWalsh. Cackett fc Bernardinil (I 2 OI 5 I I 
that these variations are not statistically significant. 

^ In m any previous studies, abundances from lAnders &: Grevessj 
1 I 1989 I I have been used for a bsorp tion models. Using these 
abundances, iBahramian et 3] ||2014| 1 find = 1.74^p Qg X 
10^^ cm~^. However in the sa me paper they show that using 
IWilms. Allen &: McCra^ ll200Clh abundances gives an equivalent 
value of 2.6±0.1 X 10^^ cm“^. 

® http://www.cosmic-lab.eu/tred/tred.php 


Terzan 5 CX9 


Hydrogen 

Helium 

Parameter 

Trial 1 

Trial 2 

Trial 3 

Njf (1022cm-2) 

3.1±0.2 

(2.6) 

2.9±0.2 

logT (K) 

6.12±0.01 

6.21±0.06 

e.ioio.oi 

NS R (km) 

(10) 

7±1 

(10) 

PL flux ^ 

1.1±0.3 

1.3±0.3 

0.9±0.3 

Xl/d.o.i 

1.11/102 

1.24/102 

1.07/102 

NHP 

0.20 

0.05 

0.29 


Table 6. Tests of hydrogen (NSATMOS) and helium (NSX) at¬ 
mosphere models for Terzan 5 CX9. Two trials use hydrogen at¬ 
mospheres: in trial 1, N^ is free, while the NS radius is fixed to 10 
km. In trial 2, N^ is fixed to the average value for Terzan 5, and 
the NS radius is free. Trial 3 is the same as trial 1, except using 
a helium atmosphere. Since we didn’t find any sign of variation 
in our initial fits, all values are tied between epochs. The neutron 
star mass is assumed to be 1.4 M0 in each case. Uncertainties are 
90% confidence, and values in parentheses are frozen. A: power- 
law fluxes are in units of 10“^^erg s“^ cm~^, in the 0.5-10 keV 
range. NHP is null hypothesis probability. 

To investigate the origin of CX12’s spectral variations 
more carefully, we perform fits letting different pairs of 
model parameters vary; kT and PL flux (Fit 1), Nir and 
PL finx (Fit 2), or PL F and PL flux (Fit 3; see Table [8} . In 
fit 1, variations in both kT and the PL flux are seen, and the 
fit is reasonable. However, the variations in kT appear to be 
anticorrelated with the variations in PL flux, which is unlike , 
for instance, the behaviour of Cen X-4 (I Cackett et al.llioid ). 
It should be noted that the signature of variation comes 
mainly from two data points (epochs 2011/4 and 2013/2). 
Performing a chi-squared test of constancy on temperature 
measured in all epochs gives a null hypothesis probability of 
0.55, providing little evidence for variation. Compared to fit 
1, fit 3 shows an increase of 14 in x^ (107 compared to 93, 
both fits have 96 degrees of freedom, giving a relative likeli¬ 
hood of 9x10“'* for fit 3 compared to fit 1 based on Akaike 
information criterion), indicating that changes in the power- 
law alone (both index and flux) are insufficient to drive the 
variations. However, fit 2, varying N_tr and PL flux, gave a 
similar fit quality to fit 1 (x^ = 92.6 compared to x^ = 93.4, 
both with 96 degrees of freedom). Thus, the variations may 
plausibly be caused by changes in the NS temperature, or 
changes in intrinsic absorption, along with changes in the 
PL flux. Since we do not have independent evidence of the 
inclination of CX12, it is quite possible that this system is 
nearly edge-on, and that we suffer varying obscuration levels 
at different times. 

We performed chi-squared and KS tests on individual 
lightcurves extracted from each observation for our target 
sources (excluding observation 13705, due to its short expo¬ 
sure). We found evidence of short (~1000 s) timescale vari¬ 
ations in Terzan 5 CX 21’s lightcurve, in observation 14475 
(Sept. 2012, Fig. [ 5 }. A KS test gives probability of constancy 
= 4.6X10“^ for this lightcurve, which is slightly higher than 
95% confidence, when accounting for the number of trials 
(80, for all sources in this study). Note that although there 
is evidence for presence of the PL in this epoch, the PL flux 
fraction is rather small (13j)^f %), suggesting that the ther- 
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mal fraction is also likely to vary. It is also possible that 
variations observed for this source are due to eclipse in the 
system, however we do not find evidence for eclipse in other 
epochs. No other source shows KS probability of constancy 
less than 0.06% (95% confidence). 


4 DISCUSSION 

4.1 Terzan 5 CX9: a NS with a helium 
atmosphere? 


Terzan 5 CX 9 shows an unusual spectral shape compared 
to other qLMXBs, with significantly lower flux at low ener¬ 
gies (between 0.1 and 1 keV), compared to similar sources. 
This could be due to a higher amount of absorption. Inves¬ 
tigating the reddening map of this cluster, we can rule out 
differential reddening. However, there could be additional 
intrinsic absorption, for instance if we are looking at the 
system edge-on. 

Another possibility is that this NS may have an at¬ 
mosphere composed of helium. Ultracompact X-ray bina¬ 
ries (with a white dwarf donor star, possessing no hydro¬ 
gen) are fairly c ommon among br ight LMXBs in globu¬ 
lar clusters (e.g. IZurek et al.l l2009l f. Helium white dwarf 
donors might thus be expected to lead to helium atmo¬ 
spheres on the NSs; there is evidence for low hydrogen 
content in X- ray bursts from s ome bright ultracomp act X- 
ray binaries JCumminel l200,fl : ICallowav et al.l l2008h . The 
X-ray spectra of helium atmospheres are similar to hy¬ 
drogen atmosphere s, but slightly harder (H Oman 1 1 19871: 
IHo fc Heinkd Im 09I1 . An observed spectrum will thus give 
different parameter values for He vs. H atmospheres, and 
several papers have considered the importance of atmo¬ 
sphere comp osition for attempts t o constrain the NS mass 
and radius j[SeryilIat_et alj_[20lj;_ [^tur^nu et al.l l2013l : 


iLattimer &: Steiner! 2014 Heinke et al. 20l4 l. For relatively 
low-count spectra such as CX9, which parameter is differ¬ 
ent from other NS qLMXBs-e.g., atmospheric composition 
(H or He), vs. Nff-may not be possible to differentiate. No 
helium atmosphere has yet been confidently identified on 
a NS, so any evidence in favor of a helium atmosphere is 
intriguing. 


4.2 Possible continuous accretion for NGC 6440 
CXI & Terzan 5 CX12 


NGC 6440 CX 1 is the brightest X-ray source in the clus¬ 
ter. It has shown multiple X-ray outburs ts over the last 
20 y ea ^rs (1998, 2001 2005 2010, 2015; lin ’t Zand et alJ 
]T99|;_ inT_Zand^t^ _ 2001 l Maxkwardt_&Swank 20051 : 


Patruno et al.ll20ld : lHoman. Poolev fc Heinkell2015l L In qui¬ 
escence, it shows a mostly thermal spectrum although there 
are signs of variations in the spectrum over the course of 
years. In the 2000 observation, we find that the spectrum 
requires a powerlaw component for a good fi t, but in later 
observ ations this component is not required. ICackett et al.l 
ll2005h argued that the changing normalization of the power- 
law component, with photon index'~2.5, was the main source 
of spectral variation. However, we suggest that the unusually 
high photon index used there led to the overestimation of 


the p owerlaw’s contribution. I Walsh. Cackett fc Bernardinil 
(l2015h also fit the spectra of CX 1 with similar models. They 
find that formally, either the powerlaw or thermal compo¬ 
nent can be the source of variation, and we agree with that 
finding. However, when we hold the NS temperature con¬ 
stant across the epochs while allowing the powerlaw nor¬ 
malization to vary (with photon index tied to a single value 
across epochs), the best-fit value of the powerlaw photon in¬ 
dex will increase to 2.7T0.5 (compared to 1.7±2 when the 
NS temperature also varies). This agrees with Walsh et al.’s 
fit results when assuming NS temperature fixed (the first fit 
in their Table 4), where they find a photon index of 3.1T0.5. 
The anomalously high photon index required for such a fit 
strongly argues that the thermal component may also be 
varying. 

Terzan 5 CX12 is well-fit by an absorbed two- 
component (thermal-hpowerlaw) model, with an Nrr column 
typical of sources in the cluster (2.6±0.1xl0^^ cm“^). CX12 
shows hints of spectral variations, concentrated at low ener¬ 
gies (< 2 keV; see Figure 2, lower right). Thus it’s unlikely 
that the powerlaw component is responsible for these varia¬ 
tions. Intrinsic absorption (e.g. an edge-on system, where an 
uneven accretion disk occasionally blocks the line of sight), 
or continuing thermal variation due to accretion are possible 
origins of these variations. 


4.3 Variability of thermal component in qLMXBs 

In our sample of 12 globular cluster qLMXBs, NGC 6440 
CXI and Terzan 5 CX12 are the only sources that show evi¬ 
dence of variations in the thermal component over timescales 
of years. In both cases, the thermal component appears to 
be the most plausible origin of the variations, but in neither 
case can thermal variations be proven. For NGC 6440 CXI, 
the best-fit thermal variation is only a 10% temperature 
change, while for Terzan 5 CX12 one observation appears to 
have a 20% drop. 

Terzan 5 CX 21 shows evidence of variability on 
timescales of hours (within observations), although we did 
not detect variations on longer timescales (between obser¬ 
vations). This variability may be due to variations in the 
(dominant) thermal component in quiescence. Alternatively, 
it may be caused by eclipses by the companion, if the system 
is at a high inclination angle. 

A crucial result of our analysis is that we have assem¬ 
bled an additional 7 qLMXBs (in addition to 4 others previ¬ 
ously reported, see Introduction) with multiple high-quality 
X-ray spectra showing little or no evidence for a powerlaw 
component; NGC 6266 CX4, CX5, CX6, and CX16, and 
NGC 6440 CXI (considering only the 2003 and 2009 obser¬ 
vations), CX2, and CX5. None of these objects show evi¬ 
dence for a powerlaw component comprising more than 12% 
of the 0.5-10 keV flux (except for the 2000 observation of 
CXI, which we exclude here). The upper limits on a pow¬ 
erlaw (assuming a photon index of 1.5) are, in four cases, 
<10% of the 0.5-10 keV flux, and <17% or <23% in two 
others. Using Chandra observations spaced over 9-12 years 
(6 for NGC 6440 CXI), we can constrain the thermal emis¬ 
sion from these NS to not have varied by more than 10% 
in any of these sources. Combined with the literature con- 
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Source 

Epoch 

logT (K) 

% Variation 

% Powerlaw fiux fraction (0.5-10 keV) 

Terzan 5 CX 9 

2003 

n 10+0.01 
^•-‘-'^-0.02 

<15 

7+7 

-6 

Powerlaw l'=(1.5} 

2009 

6.11±0.02 

<12 

<14 

Nh = 3.1 ±0.2 X 10^2 

2011-2 

6.12±0.02 

<15 

12 l« 

xS/d.o.f = 1.00/87 

2011-4 

6.11±0.02 

<12 

8 t? 


2011-9 

6.10±0.02 

- 

22 l“ 


2012-5 

o.iz_Q Q2 

<12 

<8 


2012-9 

6.11±0.02 

<12 

11 + 10 


2012-10 

6.13±0.02 

<17 

<8 


2013-2 

u.i^o_Q 008 

<13 



2013-7 

6.12±0.02 

<15 

10+“ 


2014 

C 1 OQ + O.OOS 

b.i^8_0 009 

<14 

<5 

Terzan 5 CX 12 

2003 

6.09±0.02 

00+31 

^'5_i6 

14+10 

-*-^-8 

Powerlaw r=1.8di0.3 

2009 


^^-18 

97+13 

^'-10 

Nh = 2.6 ±0.2 x 10^^ 

2011-2 


<48 

29+16 

zu_i3 

xS/d.o.f = 0.97/96 

2011-4 

o.uu_o 08 

- 

52l?8 


2011-9 

6 09'''°'^^ 
'J-^y_o.o3 

00+32 

^■^-18 

28l“ 


2012-5 

o.uy_o 02 

00+32 

27l“ 


2012-9 


<51 

04+17 

•^^-13 


2012-10 


^^-18 

29+16 

^^-13 


2013-2 

O-U^_0.02 

<38 

44ty 


2013-7 

0-UO_o 06 

<44 



2014 

_ 0 .oi 

iy+31 
^ ' -12 

29l6 

Terzan 5 CX 18 

2003 

6.09±0.02 

<17 


Powerlaw r=(1.5) 

2009 

6.09±0.02 

<17 

<13 

Nh = 2.7 ±0.2 x lO^^ 

2011-2 

6.10±0.02 

<20 

<7 

xg/d.o.f= 1.12/51 

2011-4 

6.08±0.02 

<15 

<14 


2011-9 

6.08±0.02 

<15 

7+8 

'-6 


2012-5 

6.07±0.02 

<12 

7+9 


2012-9 

o.u/ _o 03 

- 

<13 


2012-10 

o-o»y_0.o3 

<17 

<13 


2013-2 

6.08±0.01 

<12 

12±1 


2013-7 

o-U' _0.03 

<12 

10+^^ 

-‘-'^-10 


2014 

6.09±0.01 

<15 

7±1 

Terzan 5 CX 21 

2003 

O-0»^_0.03 

<15 

16l“ 

Powerlaw r=(1.5) 

2009 

o.uo_o 03 

<20 

iol“ 

Nh = 2.4 ±0.2 x lO^^ 

2011-2 

o-U' _0.03 

<23 

<12 

xS/d.o.f = 0.92/36 

2011-4 

6.03±0.03 

- 

17l“ 


2011-9 

6.08±0.02 

<26 

<10 


2012-5 

6.04±0.02 

<15 

<15 


2012-9 

o-o>^_0.04 

<15 

10+16 

i'5_ii 


2012-10 

b.UO_o 03 

<17 

11 + 14 
-*--*^-10 


2013-2 

6.04±0.01 

<12 

11+5 

ll_4 


2013-7 

6.06±0.02 

<20 

<7 


2014 

6.05±0.01 

<15 

<16 


Table 7. Spectral analyses of Terzan 5 targets with an absorbed NS atmosphere + powerlaw. All uncertainties are 90% confidence. 
Temperature variations are calculated based on the coldest measured temperature for each source. The “Epoch” values give the year and 
month, in cases where multiple observations per year were performed. Powerlaw fractions are fractions of the total unabsorbed 0.5-10 
keV fiux. All targets in Terzan 5 show significant evidence for the presence of a powerlaw component in their spectrum. In cases where 
the spectral quality allows us to constrain the powerlaw photon index (CX 12), we tied the index between observations and searched 
for the best-fit value, then froze the index to that best-fit value; otherwise we fixed the index to 1.5. The photon index used is listed in 
column 1. N/f values are in units of cm~^. 
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Figure 4. Stacked and smoothed Chandra/ACIS image of Terzan 5 in 0.3-7 keV band (total '^620 ks). Source extraction regions for 
our targets are shown with red solid circles, and regions for background regions are represented by dashed circles/annuli. For targets in 
complex regions, we extracted background from annuli. Dashed circles with slashes represent excluded regions in background extraction 
(i.e. contaminating sources). 


Terzan 5 CX12 



Fit 1 



Fit 2 



Fit 3 


Constant Nj^ & F 

Constant T/vs & F 

Constant Nh & Tjvs 


logT 

PL flux 

Nh 


PL flux 

F 

PL flux 

Epoche 

(K) 

( 10 -l 4 e]-g g-l cm“^) 

(10^2 cm-2) 

(10-1'^ 

erg s“^ cm“^) 


( 10 -l 4 ei.g g-l cm“2) 

2003 

6.09±0.02 


2±1 



2±1 

q -t -|-0.8 

'^•-^- 0.7 

7-I-6 

'-3 

2009 

6-08l“;°3 


5±2 

9 k+ 0-3 
^•^-0.2 


5±1 

2.4±0.7 

7-t4 
‘ -2 

2011-2 

6-07in3 


4±2 

9 y-(- 0.5 
' - 0.3 


4±2 

2.3±0.7 

6+3 

^-2 

2011-4 

O-'J^_0.08 


6±2 

o k+0.6 

'^•^- 0.4 


5±2 

1 q+U.D 

-^•'^- 0.7 

4±1 

2011-9 

6 09^^’°^ 

^•^=^_ 0.03 


5±2 

9 k+ 0.3 
^•^-0.2 


5±1 

2.4±0.6 

8+3 

2012-5 

o.uy_o Q2 


5±1 

2 ^- 1 - 0.3 
^•^-0.2 


5±1 

2.4±0.6 


2012-9 

6-07ini 


6±2 

9 k+ 0.5 
^•^- 0.3 


5±2 

1.8±0.8 

6^3 

2012-10 

6-08l°;°3 


5±2 

9 q+ 0.4 
^•'^- 0.3 


5±2 

2±1 

7+5 
' -2 

2013-2 

6 04^^’°^ 
^•^^_0.02 


7±2 

3.o1«;3 


6±1 

1.6±0.3 

6±1 

2013-7 

06 


8±2 

2.8l°;t 


8±2 

1.8±0.5 

8±2 

2014 

o.oi 


5±1 

2.6t°oi 


5±1 

2.2±0.4 

6±1 

xj/d.o.f 


0.97/96 



0.96/96 



1.12/96 

NHP 


0.55 



0.58 



0.20 


Table 8. Testing three possibilities to explain the spectral variations in Terzan 5 CX12. Fit 1: Fixing Njf and the powerlaw photon 
index to their best-fit value across all spectra {Nh = 2.6 x lO^^cm”^, F = 1.8), while letting the NS temperature and powerlaw flux 
vary between epochs. Fit 2: Allowing Nh a^nd powerlaw flux to vary across epochs, while fixing the NS temperature and powerlaw 
photon-index (logT = 6.08, F = 1.7). Fit 3: Only photon index and powerlaw flux are allowed to vary between epochs, with Nff and 
the NS temperature frozen (Njf = 2.5 x 10^^cm“^, logT = 6.05). powerlaw fluxes are given in 0.5-10 keV band. The NS mass and 
radius are frozen to their canonical values of 1.4 Mq and 10 km in all cases. All uncertainties are 90% confidence. NHP is null hypothesis 
probability. 
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Time (ks) 


Figure 5. Binned lightcurve of Terzan 5 CX 21 from observation 14475 (Sept. 2012). A KS test gives probability of constancy = 
4.6x10““^ for this lightcurve, which is slightly higher than 95% confidence, when accounting for the number of trials. It is possible that 
this variation may be due to eclipse. However we do not find strong evidence for eclipse in other epochs of this source. 


straints on 4 other sources d iscussed in ^1 (and in agree¬ 
ment with the recent work of I Walsh. Cackett fc Bernardinil 
[ 201 ^, this provides increasing evidence that the thermal 
X-ray spectral components of NS qLMXBs without power- 
law components to their spectra are not powered by con¬ 
tinuing accretion. (Of course, a clear detection of intrinsic 
variation of the thermal component in a source without a 
powerlaw component would disprove this hypothesis.) This 
contrasts with recent works that indicate that in some NS 
qLMXBs with powerlaw components, the powerlaw, as well 
as some of the thermal component, is power ed by continu¬ 
ing, v a riable, accretion onto th e NS surface (lCackett_et_aD 
201fll: feahram an et al.l l2014l: ICha ^bartv et al.l 20141 : 


D’Angelo et al. 20151: Wiinands et al.l ^14l). 

This evidence against continuing accretion for those 
qLMXBs without powerlaw components also supports the 
use of these qLMXB spectra to obtain constraints on the 
mass and radius (and thus on the equation of state) of NSs, 
since a lack of accretion indicates that the atmosphere will 
not be contaminated by heavy elements. 
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